Quantum Localization in Open Chaotic Systems 
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We study a quasi-Floquet state of a (5-kicked rotor with absorbing boundaries focusing on the 
nature of the dynamical localization in open quantum systems. The localization lengths £ of lossy 
quasi-Floquet states located near the absorbing boundaries decrease as they approach the boundary 
while the corresponding decay rates F are dramatically enhanced. We find the relation £ ~ p-V 2 
and explain it based upon the finite time diffusion, which can also be applied to a random unitary 
operator model. We conjecture that this idea is valid for the system exhibiting both the diffusion 
in classical dynamics and the exponential localization in quantum mechanics. 
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Quantum localization (QL) is one of the fascinating 
phenomena which cannot be expected in classical me- 
chanics. It has been observed in various physical situa- 
tions: Anderson localization (AL) in disordered systems 
[H, dynamical localization (DL) in chaotic systems Q, 
weak localization in dirty metals 0], and so on. The 
QL mainly originates from the interference among waves 
returning their initial condition. Even though classical 
statistical mechanics predicts that a particle in a random 
potential exhibits stochastic motion and thus gives rise to 
simple diffusion, the afore mentioned interference effect 
stops diffusion and localizes it within some characteristic 
length scale referred to as the localization length, £. Due 
to the interference nature of the QL, the coherence of the 
wave plays an indispensable role. 

The DL is a dynamical version of the AL in the sense 
that the distribution in momentum space stops diffusing 
and exhibits typical exponential localization. The DL 
was theoretically found in the quantum ^-kicked rotor 
(DKR) [1, 0] , a paradigm of quantum chaos and ex- 
perimentally realized by using ultracold atoms [6] . In the 
classical DKR the mean square deviation of momentum 
indefinitely increases linearly in time, while in quantum 
mechanics it follows classical evolution only for a short 
time. At the characteristic time scale, namely a break 
time it then begins to saturate so that eventually the 
quantum diffusion is completely suppressed. The formal 
equivalence between the DKR and the Anderson model 
has been proved [?|]- 

In principle every real quantum systems are coupled to 
the environment since no information can be extracted 
from completely closed systems. Thus it is a natural 
question to ask how the coupling of the quantum sys- 
tem to the environment modifies genuine quantum effects 
such as the QL. Recently the interest in open quantum 
systems has been rapidly growing in quantum chaos com- 
munity [1, 0, [1(3, 11 , [13, HH • In the semiclassical limit 
the openness introduces two major modifications onto the 
chaotic quantum systems; (i) the fractal repellers mani- 
fest themselves in quasi-eignstates Q and (ii) the density 



of states follows the so-called fractal Weyl's law [ll|. The 
severe deviation from the random matrix theory has also 
been reported in chaotic scattering when the dwell time 
of an incident particle is extremely short (shorter than 
the Ehrenfest time) [12 ] . 

Even far from the semiclassical limit the openness of 
complicated quantum systems has been an interesting 
issue, for example, the characteristics of lasin g m odes 
in chaotic microcavities 14, 15, M, [13 (see [H EU 



for review) and localization of light in random media 
[2fl 21, 22, l23l| . In this Letter we investigate the open 



DKR in quantum mechanical regime focusing on the 
characteristics of localization of lossy modes located near 
the open boundary. We found that quasi-eigenstates of 
the open DKR are separated into two kinds; one is the 
localized state whose localization length is almost equiva- 
lent to that of the DKR without the absorbing boundaries 
and the decay rate is determined from simple overlap 
argument discussed below. The other is a highly lossy 
mode located near the boundary, which is more strongly 
localized and of which decay rate is determined by con- 
sidering the finite time classical diffusion. We also show 
that all these observation and explanation is applicable 
to an open random unitary operator model. 
The Hamiltonian of the DKR is given as 



H = — + ksinx^2S(t - nT), 



(1) 



where k is the kick strength, and T is the time interval 
between successive kicks. We fix k — 14 and T = 1, i.e. 
K = kT = 14 implying the classical dynamics is fully 
chaotic. Note that the semiclassical limit implies k — > 00 
and T — > with kT kept constant, so that we do not 
consider the semiclassical regime. The time evolution 
of the open quantum DKR with absorbing boundaries 
is described by \ip{T)) = PU\^(0)), where U is a 
unitary time evolution operator for one period without 
absorption, and the operator P projects the wavefunction 
to the states satisfying \p\ < p Cl where p c represents the 
absorbing boundary. We choose p c = 1000 and set Ti = 1. 
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FIG. 2: (color online), (a) The momentum distributions of 
two typical exponentially localized QFS's with p ~ 502.86 
and 765.47, where tl > rg is satisfied, (b) As p approaches 
p c through 865.80, 961.04, and 982.88 (from the left to the 
right), the distribution becomes narrower. 
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FIG. 1: (color online), (a) The decay constants F of each QFS 
of the open DKR as a function of the average momentum p 
in log scale. The straight line represents Eq. (|3}. The inset 
shows r versus Ap = p c — p in log-log scale. The straight 
line represents Eq. Q. (b) The standard deviations o v as 
a function of p. The known theory for the DKR without 
absorbing boundary predicts a v ~ 50. 



This model has been extensively studied in the context 
of the fidelity decay or the Loschmit echo j2f| . 

The DL manifests itself via the exponentially localized 
Floquet eigenstate of U. In the open DKR the so-called 
quasi-Floquet state (QFS) can be analogously defined as 



PU\6) = e~ r e ia 



(2) 



Due to the openness of the system the Floquet eigen- 
values no longer lie at a unit circle in a complex plane 
(r > 0). Each QFS can be identified by its average mo- 
mentum, p = {(j) \p\ 4>). In one period T the survival prob- 
ability of the given QFS is given as e _2r so that the life 
time of the state, t l , is obtained by t l oc r _1 . 

In Fig. [1] we present the decay constant T and the width 
of the momentum distribution o~ p of QFSs labelled by p. 
Here two visible features are clearly observed: (i) For 
p much smaller than p Cl T increases exponentially and 
tJ p remains constant with some fluctuation, (ii) As p 
approaches p c . T more rapidly increases and a p linearly 
decreases. 

The first feature can be easily understood. For p <C p c 
the localized state has negligible influence from the ab- 
sorbing boundary so that the QFS is not so different from 
the original Floquet state of the DKR without absorbing 
boundaries as shown in Fig. EJa). This is also true even 
quantitatively since it is fairly good to estimate o~ p based 
upon the well-known relation a p ~ £ ~ D/2 ~ K 2 /4 (2(|, 



where D is the classical diffusion constant. It implies 
that the QFS can be described as <j> ~ exp(— \p — p\ /£). 
The decay constant T is then determined by consider- 
ing the overlap between the exponential tail of the lo- 
calized QFS and the absorbing region given as \p\ > p c : 
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-2r 



£ 1 fp° cx p( — 2 \p — p\ /$,)dp. It leads us to 



exp 



--A P 



(3) 



where Ap = p c ~p. Figure [TJa) clearly shows this expec- 
tation is correct. 

As p approaches p c , however, the absorbing bound- 
ary has a dramatic influence on the QFS's. They be- 
come much more lossy and even more strongly localized, 
which is unlikely because the strongly localized mode has 
smaller overlap with the absorbing region, i.e. becomes 
less lossy. It is worth mentioning that here we exploit 
the absorbing boundaries to open the system to the en- 
vironment. Usually the coupling to the environment in- 
troduces dissipation or decoherence, which destroys the 
coherence itself, let alone the localization. Therefore, one 
expects that the localization length increases. In this 
sense the absorbing boundaries are special. 

In a usual DKR there is only one important time scale, 
namely the break time rg, at which the diffusion stops 
so as to determine the localization of the Floquet state: 
<7p 2 ~( 2 ~ Dtb- In an open DKR, however, we should 
consider one more time scale, the life time, ri,. A clue 
comes from the fact that the crossover from constant to 
decreasing o~ p {p) takes place around tx ~ rg. 

The system undergoes its meaningful dynamics only 
for t < tl in the sense that the probability distribution 
no longer changes except overall decaying. If tl < tb 
is satisfied, therefore, the meaningful dynamics stops be- 
fore the break time is reached. It means that the classical 
diffusion plays a dominant role in the decay process since 
a particle disappears before the quantum suppression of 
classical diffusion takes place. In other words, the par- 
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ticle escapes when it diffusively arrives at the absorbing 
boundary. One can expect that the life time tx of a given 
QFS with p is determined simply from the duration time 
for which the particle travels from p to p c by diffusion: 
(Pc —p) 2 /D. This can be rewritten as 



TL 



D 



(4) 



which works quite well as shown in the inset of Fig. QJa). 
It is emphasized that Eq. is non-trivial because we 
are not dealing with semiclassical regime. The classical 
diffusive dynamics decides the decay constant even in the 
deep quantum regime when the loss is large enough. 

Considering the above argument the linear decrease of 
Op can also be understood. The meaningful dynamics 
stops at tl (< tb) so that a p 2 is determined not from 
Dtb but from Dt^. By using Eq. (j4|) we obtain 



(5) 



which is clearly seen in Fig. [ljb). 

In some sense very lossy modes near the absorbing 
boundary are not interesting since they decay so fast. 
They do not contribute to long time dynamics so that 
they form only broad peaks even in the scattering cross 
section. Sometimes such a mode, nevertheless, becomes 
of great importance; for example, a very lossy mode can 
play a dominant role in lasing operation, where an exter- 
nal energy input compensates the loss of the mode [HI]. 
It is also worth mentioning that the spatial shape of an 
individual mode is recently measured in the experiment 
of light in a random media, where the lossy mode strongly 



localized near the boundary is also observed [22], 

The main theme of our work is that the lossy QFS 
of the open DKR near the absorbing boundary shows 
rather stronger localization. In principle this is appli- 
cable to any system that exhibits diffusion in classical 
mechanics and exponential localization in quantum me- 
chanics. We consider one more example originating from 
one-dimensional Anderson model. Instead of Anderson's 
tight-binding Hamiltonian we exploit the so-called ran- 
dom unitary operator U = DS, where D mn — e m 6 mn 
with a random phase 9 m 2jJ , and an infinite dimensional 
matrix S is defined as 



S 



rt 
r 2 
rt 

-t 2 



V 



-t A 

-rt 
r 2 rt 
—tr r 2 
rt 

-t 2 



\ 



-t 2 
-rt 
r 2 

-tr 



(6) 



Here, S describes the scattering process via the relation 
<f>[ = Sij4>j, where (</>') represents the incoming (out- 
going) flux for a scatterer. In one-dimensional case two 
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FIG. 3: (color online) From bottom to up n = 10, 20, 30 are 
exploited for all figures, (a) and (b) are the same as Fig. HI a) 
and (b), respectively, except that a random unitary operator 
model is considered, (c) shows V versus Ax in log-log scale. 
The straight line represents Eq. Q. (d) presents a x versus F 
in log- log scale. The straight line represents Eq. ((5j). 



fluxes with the opposite direction at each site between 
two neighboring scatterers form two adjacent vector com- 
ponents, 4>2k and <j>2k+i 27]. The condition r 2 + t 2 = 1 
holds to ensure unitarity. For all calculation we fix 
t = 0.4. In order to control the width of the off-diagonal 
components, which roughly corresponds to the parame- 
ter K of a DKR, we simply multiply S: i.e. U n = DS n 
(n = 1, 2, 3, ■ ■ •). Note that the parameter if of a classi- 
cal DKR determines the maximum momentum transfer 
to a particle from each kick. The non-zero (i,j)th off- 
diagonal component of U implies there exists a non-zero 
transition probability between these two states: roughly 



speaking K ~ max(z — j) satisfying (i U j) ^ 0. One 
then expects that the larger n the bigger the correspond- 
ing effective K. We now introduce the projection opera- 
tor P and the QFS's of PU n are investigated in a similar 
way. 

In Fig. [3£a) and (b), we present the decay constant Y 
and the standard deviation a x , respectively, of an indi- 
vidual QFS of PU n for various n. Once again two visible 
features are observed. For x (= (4>\x\4>)) much smaller 
than x c {x is used here instead of p), as i increases, F 
exponentially increases and a x remains constant. As x, 
closely approach x c , however, T more rapidly increases 
and a decreases linearly. These observation are exactly 
analogous to those obtained in the open DKR (see Fig.QJ. 
As n increases, the crossover from constant to linearly de- 
creasing a x is reduced since the bigger n the larger the 
effective K, consequently the longer the break time. Fig- 
ure [3] (c) and (d) reconfirms the relations (J4j) and ([5]), 
respectively, i.e. T cx Ax~ 2 and a oc r~ 1//2 . It is empha- 
sized that these results are independent of n. 
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of quasi-bound states in light in random media. 
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FIG. 4: (color online) Husimi distribution functions of the 
QFS for (a) p ~ 502 and (b) p ~ 972. The black curve 
represents the classical unstable manifold. 



A final remark is in order. We have shown that the 
results obtained in the DKR can be directly applied to 
those of the random unitary operator model. It is noted, 
however, that chaos is different from random motion. 
First, the chaotic dynamics has more structures in phase 
space, e.g. stable and unstable manifolds associated with 
periodic orbits. In our DKR one expects they play no 
important role since we do not consider the semiclassi- 
cal limit. Indeed we deal with deep quantum regime. 
Nevertheless, the classical diffusion has a significant con- 
tribution to the decay constant when the loss is large 
enough. Even in the deep quantum regime some reminis- 
cences of unstable manifolds or chaotic repeller manifest 
themselves in Husimi plots of the lossy modes as shown 
in Fig. SJb). However, the mode far from the absorbing 
boundary does not show any similar underlying classi- 
cal structure as shown in Fig. [4ja). Note that such a 
structure does not exist in the random unitary operator 
model. Secondly, in chaotic system there is one more im- 
portant time scale, which is an Ehrenfest time te defining 
classical to quantum crossover 



12j. When r L < t e , the 



characteristics of the corresponding QFS can be strongly 
modified 12] . In our case te is extremely small so that 
one can safely ignore it [28| . 

In summary, we have shown that the QFS's with their 
life time smaller than the break time, i.e. tl < Tb, ex- 
hibit rather stronger localization and considerable loss. 
In this case, the main mechanism of the decay is deter- 
mined from classical diffusion, which gives the relation 
0. In addition, before the break time is reached the 
classical diffusion effectively stops so that the state with 
much narrower momentum distribution is formed. The 
width of the distribution is then described as Eq. ([5]). 
Such a simple explanation can also be successfully ap- 
plied to a random unitary operator model. We believe 
that our theory is valid for various physical situation 
that the diffusion takes place in classical mechanics and 
strong exponential localization exists in quantum me- 
chanics. We hope that our expectation is experimentally 
proven e.g. by direct observation of geometrical shapes 
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